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Abstract 

We study a porous medium equation with nonlocal diffusion effects given 
by an inverse fractional Laplacian operator: 

d t u - V ■ (uVp) = 0, p= (-Ay s u, < s < 1. 

We pose the problem for x £ M n and t > with bounded and compactly 
supported initial data, and prove existence of weak and bounded solutions that 
propagate with finite speed, a property that is not shared by other fractional 
diffusion models. 



1 Introduction 

We study a nonlinear diffusion model with nonlocal effects described by the system 

(1.1) d t u = V ■ (uVp), p = K,(u). 

Here, u is a function of the variables (x,t) to be thought of as a density or con- 
centration, and therefore nonnegative, while p is the pressure, which is related to u 
via a linear positive operator /C, which we assume to be the inverse of a fractional 
Laplacian. To be specific, the problem is posed for x G IR n , n > 1, and t > 0, and we 
give initial conditions 

(1.2) u(x,0) = u (x), xeW n , 
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where Uq is a nonnegative and bounded function in M. n with compact support or fast 
decay at infinity. 

The model arises from the consideration of a continuum, say, a fluid or a population, 
represented by a density distribution u(x, t) that evolves with time following a velocity 
field v(x,t), according to the continuity equation 

u t + V • (uv) = 0. 

We now assume that v derives from a potential, v = —Vp, as happens for instance 
in fluids in porous media according to Darcy's law, and in that case p denotes the 
pressure. But potential velocity fields are also found in many other instances, like 
Hele-Shaw cells. 

We still need a closure relation to relate u and p. In the case of gases in porous media, 
as modeled in the 1930's by Leibenzon and Muskat [2"5| I27j. the closure relation takes 
the form of a state law: p = f(u), where / is a nondecreasing scalar function, which 
is linear when the flow is isothermal, and a higher power of u if it is adiabatic. The 
linear relationship happens also in the simplified description of water infiltration in 
an almost horizontal soil layer according to Boussinesq. See [31] for a description 
of these and other applications. Summing up, we get the standard porous medium 
equation, u t = cA(u 2 ), or more generally, u t = A(u m ) with m > 1. 

In this paper we propose to consider the case where p es related to u through a linear 
fractional potential operator, /C = (— A)~ s with kernel K(x,y) = c\x — y|~(™~ 2s ) (i.e., 
a Riesz operator, cf. [2^1 [33] , see also Appendix for precise definitions and some com- 
ments). The interest in using fractional Laplacians in modeling diffusive processes 
has a wide literature, especially when one wants to model long-range diffusive interac- 
tion, and this interest has been activated by the recent progress in the mathematical 
theory. This literature is mostly elliptic, cf . [6J [TTJ [31] , but cf . works like |U [28] for 
related parabolic problems. 

More generally, it could be assumed that K is an operator of integral type defined 
by convolution on all of M", with the assumptions that is positive and symmetric. 
The fact the K is a homogeneous operator of degree 2s, < s < 1, will be important 
in the proofs given below. An interesting variant would be K, = (—A + cl)~ s . We are 
not exploring here such extensions. 

Extreme cases. (1) If we take in our model s = 0, so that K = the identity 
operator, we get the standard porous medium equation with m — 2, whose behavior 
is well-known, see [2j [31] for the mathematical theory and the applications. 
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(2) In the other end of the s interval, when s = 1 and we take /C = (—A) 1 , we get 



(1.3) 



u t = Vw • Vp — u 2 



Ap = u. 




XX 



V X V + c(t), 



u. In terms of v = —p x = J udx 



For c = this is the Burgers equation vt + w x = which generates shocks in finite 
time if we allow for u to have two signs. 

As a related precedent we may mention the model studied by Lions and Mas-Gallic 
[26] . who are interested in the regularization of the velocity field in the standard 
porous medium equation by means of a convolution kernel. They get a system which 
is formally like ours, but there is a big difference in the study since they assume the 
kernel to be smooth and integrable, and in fact an approximation of the Dirac delta, 
in other words, a short-range interaction. Since the kernel of the fractional operator 
(— A)~ s is k(x, y) ~ \x — y\~^ n ~ 2a \ i. e., a long-range interaction, we are far away from 
that situation, but it may serve as a previous regularization step. 

A model from superconductivity arises in recent work by Ambrosio and Serfaty 
PQ describing the evolution of the vortex-density in superconductor modeling. The 
system is similar to our system with JC = (— A) _s , s = 1 and their mathematical 
tools are quite different. On the other hand, the equation with s = 1/2 has been 
proposed by Head [21 J as the equation of motion of a dislocation continuum, and 
then u is the dislocation density and the space dimension is n = 1. The mathematical 
investigation of this case is performed by Biler et al. in [8], though in terms of the 
integrated equation v t + |f x |A(f ) = 0, where A is the Levy operator of order 1, which 
is equivalent to (— d 2 x ) 1 / 2 . 

Models of this kind arise in other contexts. Some variants will be indicated at the 
end of the paper. 

Organization of the paper. Section [2] derives the basic estimates in a formal 
way. The proof of existence of a weak solution proceeds by approximation, whereby 
the degeneracy of the equation is eliminated, diffusion is added and the kernels are 
regularized. This is technically delicate, so we first prove existence of weak solutions 
for the approximate problems posed in bounded domains in Section [3} and at the 
time basic estimates are rigorously derived. In Section H] weak solutions of the original 
problem are constructed in the whole space by passage to the limit after a tail control 
step based on a novel argument with so-called suitable "true upper barriers", cf. 
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Theorem 14.11 Such barrier method in new and turns out to be well adapted to obtain 
comparison results in the presence of nonlocal operators. 

We then establish the main properties of the solutions: Section [5] establishes the 
property of finite propagation, which is a main feature of porous media equations 
and gives rise to the appearance of a free boundary. We discuss the persistence of 
positivity in Section El 

Two sections close the paper: the Appendix, Section [7J gathers some useful def- 
initions. A final Section [8] contains comments on variants, extensions or ongoing 
work on the topic of this paper: this refers in particular to the pending questions of 
uniqueness, smoothness or asymptotic behaviour. 

Notation. We will use the notation L s = (— A) s with < s < 1 for the frac- 
tional powers of the Laplace operator defined on smooth functions in M. n by Fourier 
transform and extended in a natural way to functions in the Sobolev space H 2s (M. n ). 
Technical reasons imply that in one space dimension the restriction s < 1/2 will be 
observed. The inverse operator is denoted by fC s = (— A)~ s and can be realized by 
convolution 

K s = K s *u, K s (x) = c(n, s)\x\- (n - 2s) . 

as described in the appendix. fC s is a positive self-adjoint operator. We will write 
Hs = K l J 2 which has kernel K s ti- The subscript s will be omitted when s is fixed 
and known. For functions that depend on x and t, convolution is applied for every 
fixed t with respect to the space variables. We then use the abbreviated notation 
u{t) = «(•, t). 

2 Basic estimates 

The existence theory of weak solutions needs a lengthy process based on several 
approximations and passage to the limit that may obscure to the reader the main 
properties of the solutions. These are however very clear from usual considerations 
in mathematical physics, and the authors think that it useful for the reader to have 
them in mind as a goal. Therefore, we do at this stage formal calculations, assum- 
ing that u > satisfies the required smoothness and integrability assumptions and 
decreases fast enough as \x\ — >■ oo. The calculations are to be justified later by the 
approximation process. We fix s 6 (0, 1) and put /C = (— A) _s and H = K}l 2 . 
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Conservation of mass 



(2.1) ^- [ u(x,t)dx= [ V • {uKu)dx = 0. 

• First energy estimate: 

(2.2) f u(x, t) log u(x,t)dx = - [ (Vu • V/Cu) dx = - f \WHu\ 2 dx , 

where we use the fact that JC = H 2 , and T-L is a positive self-adjoint operator that 
commutes with the gradient. 

• Second energy estimate: 

- — j \Hu(x,t)\ 2 dx = [ T-Lu{Hu) t dx — 

(2.3) r , J^ n , 

\ ICuu t dx = / (Ku) V • (uVKu) dx — — I u\ V/Gu| 2 dx. 

JM. n JK. n JM. n 

• L°° estimate. We prove that the L°° norm does not increase in time. 

Sketch of the proof. At a point of maximum of u at time t = to, say x = 0, we have 

u t = Vu ■ VP + u AK(u). 

The first term is zero, and for the second we have — A/C = L where L = (—A) 9 with 
q = 1 — s so that 

A/Cu(0) = -Lu(0) = ~c f ^L=^dy<0, 

ji" \y\ 

where c(s, n) > 0. 

• Conservation of positivity: u > implies that u(t) > for all times. The argument 
is similar. 

• We derive next the L p estimates, 1 < p < oo : 

— / u p (x, t) dx = p j u p ~ l u t dx = —p{p — 1) / /'Vm ■ V/Gu dx = 
dt J Jr™ Ji» 

(j) - 1) J J (u(x) p - u(y) p ) AK s (x - y)u(y) dxdy = 
-(p - 1) J J {u{x) p - u{y) p ) AK s (x - y){u{x) - u{y)) dxdy 
+ (P - 1) Jf Hx) p - u(y) p ) AK s (x - y)u(x) dxdy = -I, + I 2 . 



Since u > and AK S > for < s < 1, the first integral l\ is positive. But symmetry 
means that I2 = Ii/2. We conclude that J u p dx is decreasing in time. See a related 
calculation in [15] and [T8] . 

• A standard comparison result for parabolic equations does not seem to work. This 
is one of the main technical difficulties in the study of this equation. In fact, we will 
find special situations where some comparison holds. Some partial comparison allows 
us to prove two main results of the paper. 



3 Existence I. Smooth approximate solutions 

We want to solve the initial- value problem for the equation 

(3.1) dtu = V(uV/Gi), IC = {-A)- S , 

posed in Q = W 1 x (0, 00) or at least Q T = R n x (0,T), with parameter < s < 1. 
We will take initial data uo(x) > 0, u G L 1 (R n ). We assume mostly for technical 
convenience that uq is bounded, and in the next section we will also impose decay 
conditions as \x\ — > 00. 

We want to obtain a suitable weak solution u(x,t) defined in Q. We approach 
this problem by a process that consists of regularization, elimination of the degener- 
acy, and reduction of the space domain. Once the approximate problems are solved, 
estimates are obtained that allow to pass to the limit step by step in all the approxi- 
mations to obtain in the end a weak solution of the original problem. 

Definition. We say that u is a weak solution of equation (13.11) in Qt = x (0, T) 
with initial data u G L^W 1 ) if u G L^Qt), K(u) G L^O.T : W^(R n )) } and 
u V/C(w) G L 1 (Qr), and the identity 

(3.2) J J u((f> t - V/C(n) • V0) dxdt + J u {x) <f>(x, 0) dx = 

holds for all continuous test functions in Q T such that V x is continuous, and (f) 
has compact support in the space variable and vanishes near t = T. 

[3j,l. The modifications that we use as a starting point are as follows: regularization is 
done by adding Laplacian term plus a kernel smoothing; the degeneracy is eliminated 
by raising the u = level in the diffusion coefficient. Specifically, we take small 
numbers £, 8, \i G (0, 1) and consider the equation (E(e, S, /i)): 

(3.3) ut = SAu + V -{d{u)VlC e (u)), 
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posed in Qt,r = {x E B r (0), < t < T}. A simple option for d{u) is d{u) = u + ji 
with a small fi > 0. Another option would be d{u) = \i for < u < fi and d{u) = u 
for u > fi (we prefer the former one). Besides, the equation is posed in the spatial 
domain B R = B R (0) for < t < T. We also take initial conditions 

(3.4) u(x, 0) = u (x) x G B R {0), 

where uq = uo, £ ,h is a nonnegative, smooth and bounded approximation of the initial 
data uq > 0. Finally, we take boundary data 

(3.5) u{x,t)=0 for \x\ = R,t>0. 

Let us explain now a convenient approximation of the kernel: JC £ u(t) = ( £ ~ku(t), and 
( £ is a smooth approximation of the Riesz kernel k s (x) = c\x\~( n ~ 2s ' ) corresponding 
to the inversion of the s-Laplacian on M n . In our implementation it acts on the 
extension of u(x,t) to the whole domain x G M n , which is done in the natural way, 
i. e., putting u = for |x| > R and t > 0. This approximation process is a bit similar 
to the approximation of the porous media performed by Lions and Mas-Gallic in [26J , 
but their kernel was just a mollifying kernel representing short-range effects and the 
consequences of a long-range kernel with slow decay at infinity are quite different. 

The existence and uniqueness of a solution u(x,t) = u Et s,^R(x,t) for the model is 
then more or less standard, and the solution is smooth. In the weak formulation we 
have 

(3.6) J J u(4> t - 5A4>) dxdt -JJ d(u)(VK £ (u) ■ V0) dxdt + J u (x) <j>(x, 0) dx = 

with double integrals in Qt,r = B R (0) x (0,T), and valid for every smooth <p that 
vanishes at the lateral boundary and for all large t. It is also clear that a priori 
estimates like the ones in the previous section apply to this model. In particular we 
have 

(El. a) An easy computation gives the decay of the "total mass" 



(3.7) / u(x,t)dx< / u (x) dx. 

J B R J B R 

Of course, we lose the expected mass conservation of the original problems because of 
the zero Dirichlet conditions of our present problem, but the estimate is still useful. 

(El.b) We also need conservation of nonnegativity, which is easy. 



7 



(El.c) The Lf bound is conserved, < u(x,t) < ||wo||oo, and the argument is as in 
the previous section. As a consequence, the solutions u(-,t) at time t belong to all 
L P (B R ) spaces with norm that is independent of the parameters 5, e,p and R. 

(El.d) We now introduce a version of the first energy inequality of previous section. 
We select the function F defined by the conditions F"{u) = l/d{u) and F(0) = 
F'(0) = 0. After some integrations by parts we will get 

(3.8) 4"/" F(u)dx=-5 [ ^-dx- [ \VU £ u\ 2 dx , 
K ' dtJ BR y ' J Br d(u) J Br 1 1 

1 /2 

where T-i £ = JC £ and we have used the fact that F'(u) = on E = {\x\ = R} x [0, T] 
to annihilate the boundary term in the integration by parts. This formula implies 
that for all < t < T: 

(3.9) / F(u(t))dx + 5 [ [ ^p-dxdt+ [ [ \VU £ u\ 2 dxdt = f F(u )dx. 
J b r Jo J b r d{u) J J Br J 

This implies estimates for | VH £ u \ 2 and 5\ Vu\ 2 /d{u) in L l xt (Q T)B ) and the bounds 
for such norms are independent of e, 5, R, and T. They do depend on [i > through 
the value of F = F^. Indeed, the explicit formula for F(u) is defined for u > as: 

(3.10) F» = (u + fi) log(l + {u/y)) - u, F'^u) = log(l + (u/fjt)). 

This offers difficulties when \i — y 0. (In case we take the second option, we would 
have 

Ffxiu) = {l/2jj)u 2 for < u < fi, F^iu) = u\og{u/ji) + (/u/2) otherwise, 
which is no better) 

Note. We take as JC £ the operator obtained by convolution with a standard molli- 
fication of k s of the form ( £ = k s * p £ where p £ (x) = s~ n p(x/e) and p is a C^°(lR n ), 
p > 0, p radially symmetric and decreasing; moreover, if p = a * a where a has the 
same properties. Then, we can write % £ as the operator with kernel k s / 2 * o £ . 

IS 2. We have to pass to the limit in four parameters: 5,e,p and R. The last two 
limits are the most delicate. In order to examine the convergence arguments, we can 
try to pass to the limit e — > as a next step to obtain a solution of the equation 

(3.11) u t = SAu + V ■ (d(u)V]C(u)), 
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with same initial and boundary data. Using in a precise way the above estimates we 
get convergence of u £ — > u in Lf{L^) weak. This is enough for the limit of the first 
integral of the formula of weak solution. The second integral contains the product 
u V/Gu and we need to study better the consequences of the estimates: 

(i) Since K,{u) = H{H{u)) and V/C(u) = VH{H{u)) = H(V(Hu)), we derive from 
the bound for VHu in L 2 xt the needed estimates for JC(u). Recall that H(u) and K,(u) 
are defined on all of M. n . We recall that VH(u) is a "derivative of order 1 — s of w", 
and since u is bounded, u G L™ t , we conclude that u G L 2 H^ C . By potential theory, 
it is then clear that K,{u) G L 2 H X +S . 

(ii) We also need some continuity in time. Using the equation, which expresses u t as 
the divergence of u V/Gu, together with the boundedness of u and the bound for fC{u) 
in Hl +S , we conclude that u t G L 2 (H~ 1+S ). Hence, the family of approximations to 
u is relatively compact in the sense of the parabolic compactness theorems by Aubin 
and Simon, cf. [TJ, [32] • This means that a limit u exists (along suitable subsequences) 
and u G C{[0,T] : L 2 {B R )). 

All of this is used in passing to the limit of the term jj tt(V/Cu)V0 dxdt as follows: 
we have the convergence of u £ in C([0,T] : L 2 (Br)) together with the weak conver- 
gence of p £ = K. £ u £ and Vp £ . So we can pass to the limit in this term. The conclusion 
is that we have obtained a weak solution of the initial value problem for the equa- 
tion we have mentioned, posed in Qr with zero Dirichlet boundary conditions. The 
regularity of u, Hit and Ku is as stated before. We also have the energy formula 

(3.12) f F^u{t))dx + 5 [ [ ^^-dxdt+ [ [ \VUu\ 2 dx = [ F M (uo) dx. 
Jb r Jo Jb r d(u) J J Br J 

Remark. If we pass also to the limit 5 — > 0, which is feasible with no extra effort, 
then we lose the H l estimates for u, and besides we have a problem with the boundary 
data that is maybe important. Therefore, we will keep the term 5 Au for the moment 
to avoid the problem. 

[3J.3 We will now try to pass to the limit, either as /i — > or as R — > oo, the order 
depending on convenience. In that sense we recall the detailed form of the energy 
identity 

8 ( I r dxdt + ( I I VHu\ 2 dxdt + f (u — u(t)) dx 
Jo Jb r d{u) J J Br J Br 

[ (u(t)+^)log(l + ^-)dx= I (u + /i)log(l + — )dx. 
Jb r V Jb r V 
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(Note: we write u but we should write ur since the solutions changes here with the 
radius of the ball). Terms 1, 2 and 3 are positive and behave well in both limits. It 
remains to examine the behavior of Terms 4 and 5. Both are nonnegative which is 
good for us in the case of Term 4, but Term 5 diverges as — > 0. Therefore, our 
choice is letting R — > oo. Now Term 5 is bounded by hypothesis (though the bound 
depends on /i). In the limit we easily get a solution of the problem in the whole space, 
for the equation 

(3.13) u t = SAu + V • {{u + /i)V/C(w)), xER n , t>0. 

The limit 5 — > offers then no difficulty if one wants to take it. However, the limit 
ji — > needs some extra properties. 

(3].4. One of such useful properties is the Conservation of Mass, that we establish next 
for the solutions of (13.131) . 

Lemma 3.1 Under the assumption that icq E L 1 (R n ) nL°°(IR n ) the constructed non- 
negative solutions of the previous problem satisfy 

(3.14) J u(x,t)dx = J Uq{x) dx forallt>0. 

Proof. Recall that we assume s < 1/2 if n = 1. We integrate against a cutoff 
function to get ip E C°°(M n ) supported in R < \x\ < 2R with ip = 1 for \x\ < R. We 
get 

J u t ip dx = S J u Aip dx — J (u + ^(V/Cw • Vy?) dx = I\ + 

For the typical cutoff choice we estimate the first integral as l\ = 0(R~ 2 ) using the 
fact that u(t) G L 1 (IR n ) and then I\ — > as R — > oo. As for the last integral, we do 



JCu VuV(f dx + J ulCu Alp dx + fi J JCu A<p dx . 
The latter integral can be estimated as 

/ 23 = /i f u (lCA(p)dx = //|M|iO(iT 2(1 - s) ) 



(where we use the fact that /CA has homogeneity 2(1 — s) as a differential operator) 
and this goes to zero as R — > oo. 
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Before we estimate the other two integrals we split the kernel /C into a bounded part 
Ki = min{l, ICf} and the rest K 2 — )Cf — K\ which is supported in a small ball. 
Both parts are nonnegative and moreover K\ G L°° and K 2 G L 1 . It means that 
K\ * u G L q * L p for q > qo = n/(n — 2s) (recall that we always have 2s < n) and for 
every p > 1, hence K\u G L p for all p > qo, while K 2 * u G L q * L p with q < qo and 
p > 1 which means that K 2 * u G L p for all p > 1. We conclude that /Cm G L p for 
p > q - We then have 



'22 



/c 
ulCuA(pdx< — ||ti|| g ||/CM|| 3 



which has in particular a bound of the form |/ 2 2| — C||w||i-R 2 that goes to zero as 
R — > oo. Finally, 

7 2 i = / /Cm (V« ■ Vv?) dx. 



Since Vm G L 2 and V<£> = 0(R x ) and Vv? G L p with p > n we only need /Cm G 
with a small q < 2n/(n — 2) which is true since go < 2n/(n — 2), i.e., 4s < n + 2. 

In the limit R — » oo when </? = 1, we get (I3.14p . □ 

Theorem 3.2 Let u Q G L 1 (M n ) nL°°(IR n ), u > 0. T/ien there exists a weak solution 
u — u s,e,ti of the approximate equation (13.111) posed in Qt with initial data uo, and 
u G L~(0,oo : L^W 1 )), u G L°°(Q), VH{u) G L 2 (Q). Moreover, for all t > we 
have 

(3.15) / u(x,t)dx= / u (x) dx, 



and ||w(t)||oo — 1 1 M o I loo- The first energy inequality holds, in the form 
5 [ [ ^^-dxdt+ [ [ \VUu\ 2 dxdt + 

JO JR n U + fi Jo J]&n 

/ u{t) log(w(t) + /i) dx + /i / log(l + (u/fi)) dx < 
Jm n Jr" 

/ m log(«o + /") dx + fi / log(l + {uq/ii)) dx. 
Jm n Jm.™ 

Note that the last integral is less than J uo dx. The parameters 5, e and /i are larger 
than 0. Passing to the limits 8 — > 0, e — > offers now no difficulties. But we recall 
that before taking those limits the solution is smooth in x and t, and this may be 
convenient in justifying calculations to be done below. 
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4 Tail control. Existence of weak solutions 



We start the section by stating the main existence theorem that will be proved here 
by passage to the limit fi — > in the construction of the last section. The limit that 
we want to take is not trivial because of the possible behaviour of the constructed 
solutions for large |x| which affects the convergence of some of the resulting integrals. 
We will pass to the limit in the last formula for the energy using as an extra tool a 
nice decay at infinity so that the term f u log - [u + \i) dx will be bounded uniformly 
as 11 — > 0. 



Theorem 4.1 Let u E L 1 (W l ) n L°°(M n ) be such that 

(4.1) < uq(x) < Ae~ alxl for some A, a > . 

Then there exists a weak solution u of Equation (13. ip with initial data u . Besides, 
u E C([0, oo) : L x (]R n )) ; u E L°°(Q), VU{u) E L 2 {Q). For all t > we have 

(4.2) / u(x,t)dx— / Uq(x) dx, 

and Hoe < 1 1 Mo 1 1 oo- The solution decays exponentially as \x\ — > oo as explained in 
the next tail control subsection. The first energy inequality holds in the form 

(4.3) / / \VUu\ 2 dxdt + f u(t) log(«(t)) dx < [ u log(u ) dx , 

JO JR n JR n JM. n 

while the second says that for all < t\ < t% < oo 

(4.4) f 2 [ u\V!Cu\ 2 dxdt + - [ \Hu(t 2 )\ 2 dx < - [ \H{u{t x ) | 2 dx . 

We recall that Q = R n x (0, oo), and we take < s < 1 for all n > 2, < s < 1/2 
for n — 1. The result is first proved for the solutions of the approximate equation 
with 5,e, fi > constructed in Theorem 13.21 and the information is then transferred 
to the original equation once we show that we can pass to the limit and the estimates 
are uniform. 
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4.1 Tail control 



In order to prove the above existence theorem we need to estimate a certain rate of 
decay of our solutions as \x\ — » oo. This will of course depend of a similar assumption 
that we are imposing on the initial data. In the case of the PME a possible proof pro- 
ceeds by constructing explicit weak solutions (or supersolutions) with that property 
and then using the comparison principle, that holds for that equation. Since we do 
not have such a general comparison principle here, we have to devise a comparison 
method with a suitable family of barrier functions that work because they are some 
kind of "exaggerated supersolutions" . What we need is to make sure that they do not 
admit a first contact from below. We will give to functions with such a property the 
more serious name of true supersolutions. This original technique has to be adapted 
to the peculiar form of the integral kernels involved in operator JC. We prove two 
kinds of results, the stronger one for small s. 

Theorem 4.2 Let < s < 1/2 and assume that our solution u is bounded < 
u(x, t) < L, and uq lies below a function of the form 

(4.5) U (x) = Ae~ alxl , A, a > 0. 

If A is large then there is a constant C > that depends only on (n, s, a, L, A) such 
that for any T > we will have the comparison 

(4.6) u(x, t) < Ae ct ~ alxl for all xeR n and all <t<T. 

Proof. In order to have enough regularity in the comparison argument below we 
make the proof for the solutions constructed in Theorem 13.21 in the whole space with 
parameters 6, e and /i > and we will show that the constants in the upper estimate 
are uniform with respect to such parameters if the mentioned parameters are small. 

• Reduction. By scaling we may put a = L = 1. This is done by considering instead 
of u the function u defined as 

(4.7) u(x,t) = Lu(ax,bt), b = La 2 ~ 2s , 

which satisfies the equation u t = SiAu + V • (d(u)VK{u)) with 8\ = a 2s 5/L. Note 
that then u(x, 0) < A\ e _ ' z ' with A\ = A/L. A simple calculation then shows that 
C(a, L, A) = La 2 ~ 2s C(l, 1, A/L). 
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• Contact analysis. Therefore, we assume that < u(x, 0) < 1 and also that 

u(x,0) < Ae~ r r = \x\> 0, 

where A > is a constant that will be chosen below, say larger than 2. Given 
constants C, e and r\ > we consider a radially symmetric candidate to upper barrier 
function of the form 

(4.8) U(x,t) = Ae Ct ~ r + eAe r,t , 

and we take e small. Then C will have to be determined in terms of A to satisfy a 
"true supersolution condition" which is obtained by contradiction at the first point 
(x c , t c ) G Qt of possible contact of u and U . Note that if there is contact it cannot 
happen only at \x\ = oo since u is integrable and U converges to a positive constant. 
This is the reason to add the correction term in e (use of the correction term can be 
avoided if we start the comparison argument with the solutions of the approximate 
problem posed in a ball Br(0) with zero boundary data and then pass to the limit 
R — > oo as done in the previous section). The contact does not happen at x c = 
if A is large since, putting r c = \x c \, we have at the contact point equality U = 
Ae ctc ~ Tc + eAe vtc = u and also u < 1, so that A < Ae ctc < e Tc so r c is big if A is. 
We need at least A > 1. Note that we are assuming < t c < T and T fixed (in this 
argument). 

At the point and time of contact we have u = Ae ctc ~ r ' c + eAe vtc . Assuming also 
that u is C 2 smooth, a standard argument also gives 

d r u = -Ae ctc ~ r % Au < Ae ctc ~ r % u t > AC e ctc ~ Tc + erjAe^ , 

(all of them computed at the point (x c ,t c )), and the spatial derivatives of u at x c in 
directions perpendicular to the radius are zero. Next, we put p = Ku and use the 
equation in full approximate form : 

(4.9) u t = 5Au + Vu-Vp+(u + fi)Ap, 
to get the basic inequality 

CA e ctc ~ rc + eAr] < 5 A e Ctc ~ r < - A e ctc ' rc 8^> + (A e ctc ~ Tc + eA + /x)Ap. 
Cleaning up this expression, we get at the contact point the following condition: 

(4.10) C + erie^ < 5- d~p + (1 + e e 5 + /ii)Ap , 
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where £ = r c + (r) — C)t c , \i\ = (n/A) e Tc ctc and the overline indicates that the values 
of d r p and Ap are calculated at the point of contact. 

Summing up the main ideas. In ( 14.9ft we have written the nonlinear term of the 
equation involving the fractional operator in "non- divergence form" , precisely as the 
sum of a transport term involving first derivatives on p and another term containing 
the Laplacian of p. When we evaluate those terms at the contact point with the 
barrier and arrive at ( I4.10p , the latter term is not difficult in terms of global integrals 
of u and the other two first terms in (14.91) have also simple contributions. However, 
the transport term does not have a sign and its control becomes more difficult. We 
proceed with the simple cases and leave the difficult case for the end. 

• In order to get a contradiction with inequality ( I4.10p we will estimate the values 
of d r p and Ap. For n > 1, < s < 1, we use the formula 

(4.11) t(x , t) = c J^^ dy = c j^A dy , 

with some c = c(s) > 0, which produces the singular integrals 
(a i I + \ f ( U ( x + y^) -u(x-y,t))yj 

(4.12) Px .{x,t) = Cl J 2 \y\n + 2-2s d V> 

(i = 1, • • • , n) that for s < 1/2 can also be written as 

(a i o\ t + \ f (v(x + y,t) - u(x,t))yi 

(4-13) Pa^OM) = Cl / |i n+ 2-2s d V> 



\y\ 



and finally, 



f u{x + y,t) + u{x-y,t) -2u(x,t) 
(4.14) Ap(x,t) = c 2 J | y |n+2- 2B d V- 

Here and in the sequel of the proof we denote by c, q, K, Ki different absolute positive 
constants, i.,e., constants that depend only on n and s. 

• We now estimate Ap at x c . In view of inequality ( I4.10p we need a control from 
above. Using formula (14. 14[) we see that the integral on the set \y\ > 1 is absolutely 
convergent and absolutely bounded (since u is bounded by 1). Besides, to bound 
from above the part of integral on the set \y\ < 1 we will use the fact that u(t c ) lies 
below Ae ctc ~~\ x \ with tangency at \x\ = r c ; since the numerator of the integrand is the 
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discretization of the second derivative for u at x c we use the fact that it is bounded 
above by the same expression for U to get 

u(x c + y, t c ) + u(x - y, t c ) - 2u(x, t c 



2Ae ctc - rc \ 


y\ 


2 




y\ 


n+2-2s 



d v< / ... i MJ .o-,. d v 



\r,,\n+2— 2s i 
y\<\ \V\ J\y\<l 

which is bounded since s > 0. Noting that Ae ctc ~ Tc < 1 we conclude that the term 
in (I4.10p containing Ap contributes with at most K\ to the inequality. 

• Next, we estimate the transport term involving d r p. We have 

-— f (u(x + y,t) -u(x,t))(x c -y) 
drP = c 1 J | y | W+ i- 2a d y, 

where x = x/\x\ and y = y/\y\- Since we want to get a contradiction in formula (I4.10p 
and d r p appears with a minus sign, we need to estimate this term from below. This 
integral is delicate so we split the calculation into several pieces. At this moment we 
make the further assumption s < 1/2. Then the integral for \y\ > 1/2 is bounded as: 

^ — \ i f \u(x c + y, t c ) — u(x c — y, t c )\ , f 1 

\I { \y\>m(d rP )\ < c C,l dy<c —-^dy<K 3 . 

j\y\>i/2 \y\ j y >i \y\ 

Hence, this part has the desired control. 

• The integral on the ball {\y\ < 1/2} is split again into two parts. By rotation 
we may assume that x c is directed along the first axis, x c = (r, 0, • • • ,0), where 
r = \x\. Then the integral is calculated on Oi = {y : \y\ < 1/2, y\ > 0} and on 
^2 = {y '■ \y\ < 1/2, yi < 0}. The last is easily bounded since u touches U at x c and 
lies below everywhere at time t c . Hence, 

T f ( u ( x at c ) - u(x c + y,t c ))y! f (u(x c -y,t c ) -u{x c ,t c ))y x 

-IUMP) = \ y \n + 2-2s dy = | y |n+2-2 B d V 

< f\^-y^ + f 2 ^))y U y< 1 2 ^^ d y<K 3 . 

— / L. n+2-2s Ui n+l-2s tf — a 

jQi \y\ \y\ 

where in the last inequality we have used the linear approximation for U. Again, this 
is a good control. 

• Now, the difficult part. The integral on Qi (i.e., the "half integral looking outside 
near x c ") is more delicate since the difference u(x, t c ) — u(x + y, t c ) could in principle 
drop quite abruptly even at a relatively short distance from x c and this would make the 
integral very big. However, we have a miraculous control by combining the estimate 
on the Laplacian with the good part of the estimate of d r p. 
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Lemma 4.3 With the previous assumptions and notations we have 
(4.15) - I ai (8#) + hnA^P) < -\ln 2 {d^>). 

Proof. We combine the integral of — d r p with a part of the integral for Ap as follows: 

_ r {u(x c + y,t c ) - u(x c ,t c ))yi 1 f u(x c + y, t c ) + u{x c - y, t c ) - 2u(x c , t c ) 

^ ~ I I. Ir,-l_9_9e +0 / I. lr,J_9_9» ^ 



\„.\n+2— 2s O / L.ln+2— 2s 

\y\ * jQi \y\ 

and we study carefully the integrand of Y. We have a numerator of the form 

(u(x c ,t c ) -u(x c + y,t c ))yi + ^{u{x c + y,t c ) + u{x c - y,t c ) - 2u(x c ,t c )) = 

-(j ~ yi){u{x c ,t c ) - u{x c + y,t c )) + ^{u(x c - y,t c ) - u(x c ,t c )) 

Luckily, the first term is negative, hence we conclude that estimate f)4. 15j) holds. □ 

• We can now finish the contradiction argument. All these estimates allow to conclude 
that —I(d r p) < K^. We go now back to (14.101) and conclude that the inequality implies 
that 

C + e(r] - K)e r ^~ c)tc <5 + K+^£ e r % 

This cannot happen if we choose e > 4 and C = r\ > 2K, where we keep the extra 
condition that \x and 5 must be small (at least less than 1). Then e may go to zero 
as /i — > 0. 

• Finally, we worry about the regularity of the solutions. To make the above proof 
fully rigorous, we apply the argument to the solutions of the regularized problem 
where we smooth the velocity field Vp by regularizing the kernel. We have to make 
the estimates for p r and Ap when 

p(x,t)= / K £ (y)u(x + y)dy. 



where for instance K £ (y) = K(y) for e < \y\ < 1/e, K e (y) is a parabolic cap with C 1 
fit in \y\ < e, and finally K e (y) = for \y\ < 2/e. The regularization mentioned in 
Section [3] will also do. The solutions u £ to this problem have bounded speeds, they 
are smooth and bounded with smooth and bounded p x and the previous estimates for 
p r and Ap hold uniformly in e. Passing to the limits e — > 0, the previous conclusions 
hold for any weak limit solution as constructed above, cf. equation (13.131) . The extra 
limit 5 — > offers then no difficulty □ 
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Theorem 4.4 Let now 1/2 < s < 1. Under the assumptions of the previous theorem 
the stated tail estimate works locally in time. The global statement must be replaced 
by the following: there exists an increasing function C (t) such that 

(4.16) u(x, t) < Ae c(t)t ' alxl for all xeM n and allO<t<T. 

Proof. (1) In the previous proof we had to put s < 1/2 only because of the problem 
in estimating the integral for d r p on an exterior domain, away from the contact point. 
When s > 1/2 we can estimate such integral as a convolution integral between u and 
the kernel Ki(y) = yi\y\~ n ~ 2+2s x(\y\ > !)• Now, this kernel belongs to L p (R. n ) for 
all p > n/(n + 1 — 2s), hence we only need to bound u(t) in an L q (R n ) norm with 
1 < q < n/(2s - 1) to get 

(4-17) h ■= \I{\ v \>i/2}(d^)\ < ||u(t)IUI^i|| P 

Moreover, since u G L\R n ) n L°°(IR n ) we have 

lk(0ll«<lk(*)ll! /9 ll^)lir i)/9 

We know that ||w(t)||oo < 1 and u(x,t) < Ae ct e~^ hence ||w(t)||i < cAe ct . Therefore 
the term contributes 

h < KA 1/q e ctc/q 
which allows to go back to (I4.10p and get 

(4.18) C + e{rj - K)e rc+{r >- c)tc <S + K + KA 1/q e ctc/q + ^ e r % 

The contradiction argument works as before with only one big difference. Once we put 
C = rj = KA x l q the contradiction is gotten if we restrict the time so that e ctc ^ q < 2 
which happens if 

t c < Ti = c 2 /C = c 3 K- l A- l/q = Ci A- 1/q . 
We do not play with A here, only A bigger than 2 or the like. 

(2) Once we prove estimate (14.61) for < t < T\ we can repeat the argument 
for another time interval, but now with constant A\ = Ae CTl = Ae C2 , and then 
C\ = KA l ^ q . We get a valid time interval 

T 2 = Ci A~ 1/q = c A A- 1/q e- C2,q 
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where a new factor will appear every time we repeat the iterations. In this way we 
can extend the bound up to a certain time that depends on the initial data through 
the value of A. 

(3) In order to prove that the time for which the estimate is valid goes forever we 
need to improve the L l norm of u(t) by using the fact that u < 1 together with 
u(x,t) < Ae ct e~\ x \. Summing the contributions of both upper bounds, we now get 

N*) Hi < c (log(Ae ct )) n (l + A^e- 01 ) 

so that 

h < K{\ogA + Ct) n l q (1 + A- l e~ ct ) 1/q 
and then Inequality ( 14. 18ft may be replaced by 

(4.19) C + ... < K + K(\ogA + Ct) n/q (1 + A- 1 e- Ct ) X,q + ^ e r % 

where we have dropped the terms in 5, e and fi that add no novelties or problems. 
Now we may put C = If [(log A) n/q + 2] and CT X = c 2 so that 

T 1 = c 4 (logA)-^ 

At the new starting time, we have Ai = Ae CTl = Ae C2 and then 

T 2 = c 4 (log A x )- n l* = c 4 (logv4 + c 2 )- n ^ 

and so on. Since for large k we get ~ ck~ n l q and we m ay always take q > n, the 
series ^ T k diverges, so that Estimate (14. 16j) is global in time. □ 



4.2 Proof of the existence result 



We may now pass to the limit in the weak formulation of the equation satisfied by 
u = constructed at the end of previous section. The exponential decay bound 
on the solutions, which is uniform in fi allows to improve the consequences of the 
energy inequality that is now written as 



6 J j dxdt + J Q J \ vnu v\ 2 dxdt + J u lx {t)\og + {u li {t)+ fi)dx 

< J u log(w + n) dx + n j log(l + (uo/fj)) dx + J u^(t) log~(w M (£) + /x) dx. 
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where we use the notation for the solution for clarity. Since the right-hand side is 
bounded uniformly in fi we get a uniform estimate for the family IVHwJ in L 2 (Q). 
We also have a uniform bound on Vw^ in L 2 (Q) if 5 > 0. We can therefore pass 
to the weak limit in w M — > u, and then "Hw M — > Hit and VHu^ — > VHu. The same 
happens to /Gu M and V/Cw^. The weak solution is obtained and we have 

8 J j^^-dxdt + J J \VUu\ 2 dxdt + J u(t) log + («(£)) dx 

< I u \og(u ) dx + / u(t) log~(£t(t)) dx. 



Note that the term J /ilog(l + (u //i)) dx disappears in the limit by the Dominated 
Convergence Theorem since the integrand is uniformly bounded by uq. Due to the 
decay at infinity the integral f u(t) log it (t) dx is absolutely convergent. 

The constructed solution has the same L°° bound as uq and conservation of mass 
holds by the proof of Lemma 13.11 

Passing to the limit 5 — > and e — > offers no difficulty so Theorem 14.11 is proved 
but for the second energy estimate. □ 

4.3 Second energy estimate 

Let us establish the second energy estimate for weak solutions with tail decay. We 
compute formally of the approximations where everything is justified 

-— J ip[Hu(x, t)\ 2 dx = J ipl-LuHutdx = J 'H(ip'Hu) u t dx = 
J H{ipHu) V(wV/Cw) = — j u VU{vHu) ■ V/Cw dx = 



— J uV,((pVHu) ■ WKudx — J ul-Liy^pl-Lu) ■ V/Cu dx, 
so that, putting ip = 1 — tp we have 

-— j ' p\Hu{x,t)\ 2 dx + J (fu\V ICu\ 2 dx = 

J u {ipVKu - U{ipVUu)) ■ V/Cm dx-2 J UuVp- U{u VJCu) dx, 

With the properties we have for u the right-hand side must go to zero as ip — > 1 along 
the typical cutoff sequence, at least after integration in time. 
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5 Finite propagation. Solutions with compact sup- 
port 



One of the most important features of the porous medium equation and other related 
degenerate parabolic equations is the property of finite propagation, whereby com- 
pactly supported initial data uq{x) give rise to solutions u(x, t) that have the same 
property for all positive times, i.e., the support of u(-,t) is contained in a ball B r m(0) 
for all t > and R(t) is bounded on bounded intervals < t < T. 

A possible proof in the case of the PME proceeds by constructing explicit weak 
solutions with that property (and possibly larger initial data) and then using the 
comparison principle, which holds for that equation. Since we do not have such a 
general principle here, we have to devise a comparison method with a suitable family 
of upper barriers that behave as "true (exaggerate) supersolutions" . The technique 
has been presented in whole detail in the tail analysis of the previous section and is 
here adapted to the peculiar needs of bounded support. Here is the end result. 

Theorem 5.1 Assume that u is a bounded solution, < u < L, of equation ( 13.1 ft 
with K. = (— A)~ s with 0<s<l(0<s<l/2 if n = 1), as constructed in Theorems 
14.21 and \AA[ Assume thatuo has compact support. Thenu(-,t) is compactly supported 
for all t > 0. More precisely, if < s < 1/2 and Uq is below the "parabola-like" 
function 

(5.1) U (x) = a(\x\ -b) 2 , 

for some a, b > 0, with support in the ball B b (0), then there is C large enough, such 
that 

(5.2) u{x,t) < a(Ct- (\x\ -b)) 2 

Actually, we can take C = C (n, s)L( 1 ' 2 } +s aP~' 2 ^~~ s . For 1/2 < s < 2 a similar 
conclusion is true, but now C is an increasing function of t and we do not obtain a 
scaling formula for its dependence of L and a. 

Proof. The application of the method is very similar to the case worked out the 
tail control section. Therefore, we will dispense with some of the technicalities of 
regularization to gain space and clarity. We assume that our solution u(x,t) > has 
bounded initial data Wo(^) = u(x,t ) < L and also that u is below the parabola 
Uq(x) = a(\x\ — b) 2 , a, b > 0; moreover, the support of uq in the ball of radius b and 
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the graphs of uq and Uq are strictly separated in that ball. We take as comparison 
function U(x,t) = a(Ct — (\x\ — b)) 2 and argue at the first point and time where 
u(x,t) touches U from below. The fact that such a first contact point happens for 
t > and does not happen at x = oo is justified by regularization as before. We put 
r = \x\. 

By scaling we may put a — L — 1. See detail of the reduction step below. We 
examine in detail the situation in which the touching point (x c , t c ) is not the minimum, 
say, x c lies at a distance from the front := b + Ct, so that b + Ct c — \x c \ = h > 0. 

Note that since u < 1 we must have \h\ < 1. Assuming also that u is C 2 smooth, a 
standard argument gives 

u = h 2 , u T = —2h, Au < 2n, u t > 2Ch, 

all of them computed at the point (x,t) = (x c ,t c ). Putting p = Ku and using the 
equation u t = Vw ■ Vp + uAp, we get the inequality 

(5.3) 2Ch < -2/iJv + h 2 Ap~, 

where the overline indicates that the values of p r and Ap are calculated at the point of 
contact. Moreover, u(x, t c ) < {xf — x) 2 for all x e M n . In order to get a contradiction 
we will estimate the values of p r and Ap. The formulas for p, p r and Ap in terms of 
u are given in fTCTTj) . (HT2"|1 . (03|1 . 

• Estimating Ap at the contact point offers no novelties. As before and in view 
of inequality (15. 3 p we need a control from above. Using formula (I4.14p we see that 
the integral on \y\ > 1 is absolutely convergent and bounded (since u is bounded). 
Besides, to bound the integral for \y\ < 1 from above we will use the fact that u lies 
below a parabola, hence 

f u(x + y, t c ) + u(x - y, t c ) - 2u(x, t c ) f 2\y\ 2 

I , L,|n+2-2a ay ~ /, , U.|n+2-2s 

J|)/|<i 1 2/ 1 ^|y|<i 1 2/ 1 

which is bounded since s > 0. We conclude that the term in (15. 3p containing p xx 
contributes with at most Kh 2 to the inequality, where K > is an absolute constant. 

• Next, we estimate p r . Much of the argument is similar to the tail analysis but 
the end is more delicate. Since we want to get a contradiction in formula (I5.3P and 
the coefficient u r = —2h of p r is negative, we need to estimate this term from below. 
Using formula ( 14. 12ft we see again that the integral for \y\ > 1/2 is absolutely and 
uniformly bounded by a constant K 2 . 
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• The integral on the ball {\y\ < 1/2} is split into several parts. We will drop for 
convenience of writing the dependence on t c in the formulas the follow. By rotation 
we may assume that x c is directed along the first axis, x c = (r, 0, ■■• ,0), where 
r = \x\. Then the integral is calculated on Qi = {y : \y\ < 1/2, y\ > 0} and on 
^2 = {y '■ \y\ < l/2,yi < 0}. The last is easily bounded since u touches U at x c and 
lies below everywhere at time t c . As in the tail analysis we have, 

t n*-\ f ( u ( x c) -u(x c + y))y x ^ t, 
-InMp) = ^ZTs dy<K 3 . 

• The integral on Qi (i. e., the "half integral looking outside near x c " ) is more delicate 
as we have said, since the difference u(x, t c ) — u(x + y, t c ) could in principle drop quite 
abruptly even at a relatively short distance from x c and this would make the integral 
very big. There is a part with \y\ between Oh and 1 for any given 9 G (0, 1) that is 
easy (note that < u(x, 0) — u(x +y,t) < u(x, 0) = h 2 ): 



Jeh \y\ Jeh \y\ 

and this is good even for small h. 

• The last part of the integral for p r , over the half-ball H 1 = {0 < \y\ < 9h} with 
yi > 0, is in principle bad since u(x c + y) could drop abruptly, thus making the 
integral very negative or even divergent near y = 0. We are going to combine the 
integral of — d r p with a part of the integral for Ap as follows: 

_ f (u(x c + y) - u(x c )yi t h f u(x c + y) + u(x c - y) - 2u(x c ) 



U,|n+2-2s O / L,|n+2-2s 

jhi \y\ A jhx \y\ 

and we study carefully the integrand of Y . We have a numerator of the form 

(u(x c ) - u(x c + y)yi + ^(u(x c + y) + u(x c - y) - 2u(x c )) = 

h h 
-(j ~ yi)( u ( x c) - u ( x c + y)) + -^{u(x c -y)- u(x c )) 

Luckily, the first term is negative (we take 0<#<l/2(a security factor), hence we 
conclude that 

-I Hl <&p) + ^ Hi (Ap) < -2 7 n 2 (^P)- 

• We may now sum up all the terms in the right-hand side of (15. 3p and show that 
they are bounded above by Kh where K is a uniform constant. Therefore, for large 
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C inequality ( 15. 3 p is impossible, hence there cannot be a contact point with h 7^ 0. 
In this way we get a minimal constant C = Co(n, s) for which such contact does not 
take place. 

• Reduction Step. We use it to get the dependence on L and a. Here, it goes as 
follows: since the equation scaling is 

u(x,t) = Au(Bx,Tt) 

with parameters A, B, T > such that T = AB 2 ~ 2s , if we have done the proof for u 
that has height 1 and is below (\x\ — bo) 2 initially, and get a comparison with a U as 
above with speed Co > 0, then the assumptions u(x, t) < L and u(x, t ) < a(\x\ — b) 2 
initially, are satisfied if we put 

A = L, AB 2 = a, b = b /B, 

i.e., A = L, B = (a/L) 1 / 2 , and then T = a}~ s L s . The new speed is then 

C = = Coa l / 2 - s L l / 2+s . 
B 

• We still have to consider the modification of the proof when 1/2 < s < 1. The only 
problem is the estimate of d r p on the exterior of a ball. This is done as in the tail 
control case. 

We next lemmas complete the details of the comparison proof that have been left 
out in the previous lemma. 

Lemma 5.2 Under the assumptions of Prop. 15.11 there is no contact either atu = 0, 
in the sense that strict separation of u and U holds for allt>0 if C is large enough 
as in the previous lemma. 

Proof. Precisely, what we want is to eliminate the possible contact of the supports 
at the lower part of the parabola. Instead of doing this by analyzing the possible 
contact point, we proceed by a change in the test function that we replace by 

U e = (Ct- (\x\ - b)) 2 + e(l + Dt) for \x\<b + Ct, 

and U £ = e(l + Dt) for lx\ > b + Ct. Here e > is a small constant and D > will 
be suitably chosen. Assume that the solution starts at t = and touches this test 
function U E for the first time at the time t = t c . 
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The argument for a contact at the points \x\ < b + Ct c where U £ is a parabola works 
like previously. 

The argument at contact points |x| > b + Ct c leads to 
(5.4) De < e(l + Dt c )Ap~, 

where the overline indicates as before that the value Ap are calculated at the point 
of contact. Moreover, u(x,t c ) < Ui(x,t c ) for all 16K. If we are able to prove again 
that Ap < K we will get 

De < e(l + Dt c )K 

which is contradictory if D > K, say D = 2K, and t c < l/D = 1/(2K), This 
estimate is uniform in e and gives in the limit and upper bound of the form u(x, t) < 
(Ct — (\x\ — b)) 2 for all < t < 1/2K, and as a corollary, the support of u bounded 
on the right by the line \x\ = Ct + b in that time interval. After this time the process 
can be repeated and the conclusion is true for all times. □ 

• We now reflect for a moment on the regularity requirements. Arguing as in the tail 
control case by using the smooth solutions of the approximate equations, the previous 
conclusions hold for any weak limit solution. 

5.1 Consequences. Growth estimates of the support 

The following analysis is done for s < 1/2 and concerns bounded solutions with 
compactly supported initial data. By free boundary we mean, as usual, the topological 
boundary of the support of the solution S(u) := {(x,t) : u > 0}. 

Corollary 5.3 Let uq be bounded above by L with uq(x) = for \x\ < R. Then, we 
get an estimate for the free boundary points of the form \x(t)\ < R + C2 1 1 ^ 2-2 ^ if 
s < 1/2. 

Proof: We know that the support of u(-, t) is bounded for all times, say, it is contained 
in a ball of radius r(t). We take a time t\ and find a parabolic barrier as before, with 
coefficient a > that is initially above and separated from u(-,ti). This can be done 
by choosing first r\ such that ar\ = L and then putting in the formula of the parabolic 
barrier b = r(ti) + r% + e, and then we can go forever in time in the comparison. Using 
the speed estimate in Theorem 15.11 we get in the limit e — > 0: 

r(t) - r(h) - n < C(t - h) = CLihf'^a^it - h) = C L{h)(t - h)/r\- 2s 
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Here r\ > is free by moving a. We can use the L°° bound L(t\) < L(0) We get 

r(t) <r(t 1 )+r 1 + C L{t-t 1 )/r{- 2 ' 

Now take t\ = and optimize the right-hand expression in r\ > for s < 1/2. . 
Notice that in the limit s = 1/2 we would get linear growth, while for s = we get 
the standard t 1 / 2 growth of the Porous Medium Equations under these assumptions. 

6 Persistence of positivity 

We establish another property that plays an interesting role in the theory of porous 
medium equations to avoid degeneracy points for the solutions. It is called persistence 
of positivity. For continuous solutions it implies non-shrinking of the support. 

Theorem 6.1 Let u be a weak solution as constructed in Theorem (14.1 ft and assume 
that uq(x) is positive in a neighborhood of a point xq. Then u(xq, t) is positive for all 
times t > 0. 

Proof. This issue allows for another use of the technique presented in the tail analysis, 
but this time with a true subsolutions. We assume that u (x) > c > in a ball B R (x ). 
By translation and scaling we may assume that xq = and c — R — 1. The idea is 
to study the contact point with a parabola that shrinks quickly in time, like 

U(x,t)=e- at F(\x\), 

with F suitably chosen and a > large enough. Firstly, we choose F to be radially 
symmetric and decreasing, with F(0) = 1/2 and -F(|x|) = for \x\ > 1/2. The 
contact point (x c , t c ) is sought in _E>i/ 2 (0) x (0, oo). By approximation we may assume 
that u is positive everywhere so no contact at the parabolic border is assumed. At a 
positive contact point we have u t <U t = —all, Vu = e _at F / (|x|)e r and the standard 
arguments on the equation imply 

-aF(\x\)>F'(\x\)p- + F(\x\)A^ 

Now we have 

p{x,t) = Ku{x,t) = e- at {lCF){r) 
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As in the tall control analysis we prove that A x p is bounded uniformly. The novelty 
is that F' < implies that p r < so that the term F'(\x\)p^ > 0. In this way the 
inequality implies 

a < Ke~ at 

which is false if a > K. Hence, under this size assumption there can be no contact 
point, and U is a true subsolution, and positive for all times in Bi/ 2 (0). □ 



7 Appendix. Fractional Laplacians 

We collect some data on fractional Laplacians for the reader's convenience 

Fractional Laplacians and potentials. According to Stein [33], Chapter V, the 
definition of (— A)^/ 2 is done by means of Fourier series 

((_A)^/)» = {2n\xff{x) 

and can be used for positive and negative values of 0. For /3 = —a negative, with 
< a < n, we have the equivalence with the Riesz potentials 

- a /2 f _ T rt\ ._ 1 f f(v) 



7(a) Jw 1 F - V\ 

(acting on functions of the class S for instance) with precise constant 

7 («) = vr n/2 2 a r(«/2)/r((n - a)/2). 

Note that 7 — > 00 as a — > n, but j/(n — a) converges to a nonzero constant, 
7r™/ 2 2 ri ~ 1 r(n/2). Stein also mentions the Bessel potentials which are associated to 
the modified inverse Laplace operators 

{{I-A)-«l 2 f)=l a {f) 

The Bessel potential has a kernel G a (x) that is better behaved at infinity, though not 
given by a simple kernel (see [33], page 132). 



Fractional Laplacians via extensions 

In the case s = 1/2 it is the well-known technique of harmonic extension to the upper 
plane of the space with one more dimension and then taking the boundary normal 
derivative, and has been used in the study or variational problems with thin obstacles 
as in [3], see also [HIED]- For s 7^ 1/2 the method has been recently developed in [T2"] . 
it involves elliptic equations with weights and weighted normal derivatives. 



27 



8 Comments and extensions 



Extensions. Very different versions to the evolution process are obtained when the 
pressure is related to u in other ways. Thus, we can use a pressure-density relation of 
the form p = f(u) with / an increasing function, and then the model equation would 
be 

(8.1) d t u = V -{uVK{f{u))). 

Many of the results proved here should apply to this model. Another possibility 
consists of equations of the form d t u = V ■ (/(u)V/C(ii)). 

Relaxation. Chemotaxis models. We could also relax the relation of p to u into 
the form 

(8.2) d t p + {-A) s p = u. 

This reflection is motivated by a very important system, the Keller-Segel chemotaxis 
model, [231 122] , in which the phenomenon which is modeled by u is not diffusion but 
the concentration of a certain population and p is replaced by variable c proportional 
to the concentration of the chemical substance responsible for the aggregation of the 
population. A suitable general system is proposed in the form 

(8.3) u t = sAu- V • (uVc), Set + K' l c = f{u, c). 

The standard chemotaxis model uses K = —A and f(u, c) = u — be. In the limit case 
where 5 and b are zero, and if we use as K, an integral operator as described above, we 
get a model with a term like ours but note the different sign, u t = eAu — V- (uVJC(u)), 
which is a consequence of the fact that we are dealing with aggregation and not 
diffusion. The study of these equations is also of interest. 

Finite propagation. The finite propagation property is not true for other alterna- 
tive models of porous medium equation with fractional diffusion like the ones studied 
in [1], [28] The last reference deals with the model 

(8.4) _ + (_A) 1 /2 ( | ur i M) = o, 

For any m > m* = (n — l)/n, it is proved that a unique nonnegative strong solution 
of this problem exists for data in Li (R n ) and is strictly positive. The maximum 
principle applies to this problem. 
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Uniqueness and comparison. These are widely open issues. We have used in 
the present paper comparison with what we call "true supersolutions" and "true 
subsolutions" . In one space dimension, a uniqueness proof has been obtained in [8] 
by integrating the equation (with respect to x) and using then solutions in the sense 
of viscosity. Such trick is not available in several dimensions. 

Evolution and regularity of free boundaries. This is quite important topic 
motivated by the property of finite propagation. 

Smoothing L 1 into L°° and C a regularity. These topics will be treated in [13] . 

Asymptotic behaviour. This topic is under study. Let us outline the main details. 
There exists a family of self-similar solutions for this problem, in the spirit of the fun- 
damental solution of the linear problems or the Barenblatt solutions of the standard 
Porous Medium Equation. The spatial profile of such solutions is obtained by solving 
for the pressure p an obstacle problem with a truncated paraboloid as obstacle; the 
corresponding density u is then the mass of the negative fractional Laplacian of p, 
supported on the contact set; all this fits perfectly into the elliptic theory described 
[6]). The details of the construction, as well as the convergence of a typical solution 
to such asymptotic profiles after suitable scaling, will be established in an upcoming 
publication, [T4"] . 
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